Structure of the Wake of a Magnetic Obstacle 
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We use a combination of numerical simulations and experiments to elucidate the structure of the 
flow of an electrically conducting fluid past a localized magnetic field, called magnetic obstacle. 
We demonstrate that the stationary flow pattern is considerably more complex than in the wake 
behind an ordinary body. The steady flow is shown to undergo two bifurcations (rather than one) 
and to involve up to six (rather than just two) vortices. We find that the first bifurcation leads to 
the formation of a pair of vortices within the region of magnetic field that we call inner magnetic 
vortices, whereas a second bifurcation gives rise to a pair of attached vortices that are linked to the 
inner vortices by connecting vortices. 

PACS numbers: 41.20.Cv, 47.65.-d, 47.80.Cb, 47.90.+a 



When a liquid metal moves relative to a localized mag- 
netic field - a magnetic obstacle - the induced eddy 
currents produce a Lorentz force that creates vorticity. 
Whereas the flow pattern around a mechanical obstacle, 
such as the cylinder shown in Fig.[TJ are well documented 
[lj], the structure of the wake of a magnetic obstacle is 
poorly understood even in the seemingly simple steady 
state. In the present Letter we demonstrate that the first 
bifurcation in the flow past a magnetic obstacle leads to 
the formation of a vortex-pair inside the magnetic field 
and that the formation of attached vortices, the analog 
to Fig. []J, is the result of a second bifurcation giving rise 
to a remarkably complex six- vortex pattern. 

Understanding the wake of a magnetic obstacle, whose 
prototype is shown in Fig|2l is of considerable funda- 
mental and practical interest. On the fundamental side 
@- H, 0, Q, a magnetic obstacle represents a system 
with a rich variety of dynamical states whose behavior 
is governed by two parameters, the Reynolds number 
Re = uoH/ v (where uo is velocity, H is the characteristic 
linear size, and v is the kinematic viscosity of the fluid) 



and the magnetic interaction parameter N = ctHBq / puo 
(where a is the electric conductivity of the fluid, p its 
density and Bo the magnitude of the magnetic field) Q . 
On the practical side, spatially localized magnetic fields 
enjoy a variety of industrial applications p exemplified 
by electromagnetic stirring, electromagnetic brake, and 
non-contact flow measurement 

An analogy between the hydrodynamic obstacle and a 
braking localized magnetic field has been supposed from 
the beginning of Magnetohydrodynamics. However, up 
to now the corresponding reverse flow behind a magnetic 
obstacle has never been observed experimentally. The 
earliest 2D calculations [8[ showed a kind of recircula- 
tion, but the especially designed experiments [9] failed to 
confirm it. Recent 2D numerical papers of Cuevas et at. 
revived the term 'magnetic obstacle' [15] and found a 
Karman vortex street [10[ for large Re, as well as a vor- 
tex dipole in a creeping flow [111 ]. 

We shall derive logically, for the first time, new possi- 
ble stationary magnetohydrodynamic flow patterns and 
illustrate them by both physical experiments and 3D nu- 
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FIG. 1: Streamlines of the steady flow around a solid cylin- 
der and corresponding streamwise velocity along the center- 
line: (a) without vortices and (b) with attached vortices. No- 
tice, the region of attached vortices corresponds to negative 
streamwise velocity. 
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FIG. 2: (color online) Magnetic obstacle (inset) and stream- 
wise velocity along the center line from 3D numerics (a — c) 
and experiment (d): N=4, Re=100 (a), N=11.25, Re=100 
(b), N=11.25, Re=400 (c), N=11.25, Re=2000 (d). Dashed 
vertical lines mark borders of the magnetic obstacle. 
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FIG. 3: (color online). Flow streamlines in the middle plane corresponding to profiles of Fig. [2] involving inner (red b — d), 
connecting (green c, d), and attached (magenta c,d) vortices. Dashed lines mark borders of the magnetic obstacle. 

merical simulations. Technical details to reproduce our 
quantitative results are given to the end of the paper. 

Let us recall that interaction parameter N represents 
the ratio of the Lorentz forces to the inertial forces and 
Reynolds number Re is the ratio of inertial and viscous 
forces. From this we derive that (i) the higher N the 
stronger is the braking Lorentz force inside the obsta- 
cle and (ii) the larger Re the more pronounced is the 
stagnant region behind the obstacle. Now we analyze all 
possible dependencies of streamwise velocity along the 
centerline (centerline curves), Fig. [2j analogously as it 
is shown in Fig. dJ Figure [3] illustrates these centerline 
curves by flow patterns obtained from 3D numerical sim- 
ulation and physical experiment. It is obvious to con- 
clude that streamwise velocity inside the magnetic ob- 
stacle is positive for small N (curve a) and might be neg- 
ative when N is higher than a critical value N™ (curves 
b — d). Past the obstacle the flow structure is analogous 
to that of a solid cylinder shown in Fig. [TJ at low Re the 
velocity is positive (curves a, 6), and it might be negative 
when Re is larger Re c (N) (curves c, d). 

It is a signal for a reverse flow when a centerline curve 
becomes negative. Therefore, we deduce that in the most 
developed case there must be a pair of vortices inside the 
obstacle, which we call inner magnetic vortices and show 
by red color in Fig.[3](6 — d), and a pair of hydrodynamical 
attached vortices past the obstacle marked by magenta 
color in Fig. [3^c, d). 

The inner vortices and the attached vortices show the 
same direction of rotation determined by the main flow 
movement. This means that a continuous flow only can 
exist, if the pairs of inner and attached vortices are con- 




FIG. 4: (color online) Topology of the flow involving inner 
(red), connecting (green), and attached (magenta) vortices, 
bold dashed lines mark borders of magnetic obstacle. Cross- 
marked are fixpoints, black lines are heteroclinic orbits. 



nected by contra-rotating connecting vortices shown by 
green color in Fig.[3tc, d). These connecting vortices cor- 
respond to a local maximum of centerline curves in Fig. [2] 
at x « 2.5. The resulting flow structure of Fig. [3][c, d) 
becomes clear from Fig. [4j where the 2D flow structure 
in the middle plane is organized around six hyperbolic 
fixpoints and their corresponding heteroclinic invariant 
manifolds. As can be derived from nonlinear topological 
dynamics [12] no other two-dimensional volume conserv- 
ing and continuous differentiate flow structure is possi- 
ble, in which two pairs of vortices with the same sense 
of rotation can coexist. The inner and connecting vor- 
tices with the surrounding invariant manifolds replace the 
rigid body of Fig. [TJ Notice, that a virtual (magnetic) 
obstacle induces real vortices inside. This is a direct rep- 
resentation of the hydrodynamic potential theory where 
a real physical body is modelled by means of a virtual 
vortex dipole flow. 
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Recirculation inside the obstacle is dependent on the 
interaction parameter N being inverse proportional to 
velocity. Therefore, starting from a stable flow without 
vortices, inner vortices and then the six- vortex pattern 
involving attached vortices appear by decreasing the flow 
rate. This is a paradoxical behavior being in contrast 
to the ordinary hydrodynamics where attached vortices 
always appear at increasing flow rate. 

To confirm the foresaid theoretical finding experiments 
have been performed. The measured centerline velocity, 
curve d in Fig. [2j agrees very well with simulated re- 
sults in the region of the magnets. Behind the magnets 
there is a discrepancy which is obviously explained by 
the fact that the experimental Re = 2000 is essentially 
larger than in 3D simulation. However, the experimental 
flow structure shown in Fig. 09(d) is qualitatively the same 
as in Fig. [31(c). Therefore, the found six- vortex pattern 
is stable in the experiments even at Re = 2000. Both 
experiment and numerical simulation show that with in- 
creasing Re and constant N the attached vortices grow 
while the connecting vortices slightly shrink. 

Direct quantitative comparison between experiment 
and simulation is not possible due to the following rea- 
sons. First, inlet flow is laminar in simulations and tur- 
bulent in experiments. Second, high experimental Hart- 
mann number is difficult to handle in a proper 3D simu- 
lation since Hartmann layers must be finely resolved on 
the bottom and top walls. For qualitative comparison 
these factors are not important, since Re number plays 
no role for the core flow [13[ and turbulent pulsations are 
suppressed by magnetic field [14J. The slight spanwise 
asymmetry of the experimental flow, Fig. [3]i, is due to a 
bend of the flow circuit in front of the channel inlet. 

Let us now shortly describe some technical issues. The 
governing equations which were solved in the present 3D 
simulation to obtain Fig. [2j [3] are composed of Navier- 
Stokes plus Maxwell equations for moving medium and 
Ohm's law. Induced magnetic field is assumed to be in- 
finitely small compared to an external magnetic field 0] . 
Then, these equations are the following: 

+ (u • V)u = -Vp + -Au + N(j x B), 
ot Re 
j = -V(/> + uxB, V-u = 0, A(/> = V • (u x B), 

where u is velocity, B is external magnetic field, j is elec- 
tric current density, p is pressure, (j) is electric potential. 
The above equations were solved by an own self-made 
3D solver for the same rectangular duct as in the exper- 
iments. Technical numerical details are given in [l3| ]. 

The experimental data presented in Fig. [2j [3] were 
measured by one of the authors, O.A., with the Ul- 
trasonic Doppler Velocimeter DOP2000, on the exper- 
imental setup described in details in [14]. This setup 
is represented by the Ga-In-Sn eutectic alloy flowing 
in the rectangular duct with sizes (Length x Width x 
Height) =2 (L x x L y x if) =(50 x 10 x 2). The mag- 
nets assembled in a yoke have horizontal sizes (Length x 
Width) =2 (M x x M y )=(3 x 4) and distance between the 



poles h = 3. All numbers are in centimeters. Magnetic 
field strength measured in the center of the magnetic gap 
is B = 0.4T. 

Re and N through the paper are defined with the pa- 
rameters H, uo and Bo equal to half- height of the duct, 
the mean flow velocity, and field intensity measured in 
the center of the magnetic gap, correspondingly. 

External magnetic field was modelled as a field from 
two permanent magnets occupying a space Q — \\x\ < 
M x , \y\ < M y , \z\ >= h}, then B(r) = V / n B d (r, r')dr', 
where Bd(r, r') = d z r(l/\r — r'|) is a field created by a 
single magnetic dipole. After algebraic calculations: 

B 7 (x,y,z) = ^2 ijkF 7 (x - iM x ,y - jM y ,z - kh), 

i,j,k=±l 

where 7 = x,y,z are for the magnetic field compo- 



nents, F x (x,y,z) = tanh 



F y (x,y,z 



F z (x,y,z) = tan 



-1 zr_ 
xy ' 



and r = (x 2 + y 2 



tanh- 1 
z 2 )- 1 ' 2 . Such 



an analytic approach is rather precise, as was tested by 
comparing the calculated magnetic field with the field 
measured experimentally. 

Let us now turn to final results presented as stability 
diagrams in Fig. They summarize the flow behav- 
ior around the magnetic obstacle dependent on the three 
parameters Re, N and ft, where ft = M y /L y is the con- 
strainment ratio between the lateral size of the magnet, 
2M y , and the width of the channel, 2L y . The foresaid 
simulated and experimental results correspond to a con- 
strainment ratio of ft = 0.4. 

In Fig 5(a)| the constrainment factor is varied for the 
Reynolds number Re = 100. The less is ft, the less is the 
influence of side walls. The case of ft — > corresponds to 
a free flow. The larger is ft, the more spanwise uniform 
is the braking Lorentz force. Therefore, in order to in- 
duce inner vortices at larger ft it is necessary to apply a 
larger critical interaction parameter N™ as one can see 
in Fig. 5(a)| At ft < 0.5 the critical value of N™ is of 
the same order of magnitude, N™ w 6, and after ft « 0.8 
rapidly increases up to N™ ^110atft = l. We believe 
that it is impossible to provoke magnetic vortices at high 
n corresponding to spanwise uniform magnetic field. 

At small constrainment ratio (kC 1) and fixed M x the 
magnetic obstacle becomes very slim (magnetic blade) 
compared to the channel width. Then, attached vortices 
are not formed because of well streamlined shape of the 
blade, while the magnetic vortices are originated as be- 
fore. Now, these magnetic vortices are elongated along 
side walls of the magnetic blade and are expanded in 
spanwise direction if interaction parameter N increases. 

Finally we present in Fig. |5(b)] the stability diagram for 
ft = 0.4. It shows, derived from a series of 3D simulation 
and experiments, bifurcation lines N™(Re) and N^(Re) 
separating the existence regions of flow without vortices 
(I) , with inner vortices (II) , and flow including inner mag- 
netic vortices, connecting vortices and attached vortices 
(III) . To depict the borders between regions I-III we have 
calculated several foothold points given in Fig. |5(b)| by 
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FIG. 5: (color online), (a): Stability of inner magnetic vortices in dependence on constrainment k and interaction parameter 
N at fixed Re = 100. (b): Global stability diagram giving the existence regions of all three flow types described in this Letter 
in dependence on (Re, N) at fixed k = 0.4. 



filled symbols. The open triangle symbol at Re = 4700 
was obtained experimentally. Bifurcation points were de- 
rived by monitoring the first (for inner) respective second 
(for attached vortices) minimum on centerline curves as 
shown in Fig. [2j These minima must be exactly equal to 
zero at the bifurcation between two flow patterns. The 
first bifurcation takes place slightly above N w 5 and is 
almost independent of Re. If the magnetic vortices are 
present, a minimum Reynolds number, Re™, is necessary 
for the second bifurcation to the appearance of attached 
vortices analogously to ordinary hydrodynamics. Since 
we have never observed attached vortices in the simula- 
tion with Re — 100 and N < 36, Re™ is given near 100. 
A series of experiments at various Re resulted in a critical 
value of Re c ~ 4700 for inner vortices which corresponds 
to N™ « 5 fitting well to the numerical values. 



In summary, new phenomena are discovered for con- 
ducting liquid flowing through a magnetic obstacle. One 
is the appearance of stable inner vortices inside the 
obstacle at magnetic interaction parameter N > N™. 
The other is a new six- vortex flow pattern appearing at 
N > N^(Re). This pattern is composed of three pairs 
of vortices: inner, connecting, and attached. All the vor- 
tices are found both in 3D numerical simulation and in 
physical experiments. 
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